McGill-96/38 
UA/NPPS-7-96 
|hep-ph/ 961 1401 

September 1996 



Direct photon production in polarized hadron 
collisions at collider and fixed target energies 

A. P. Contogouris a,fe ' 1 and Z. Merebashvili a '*' 2 

a. Department of Physics, McGill University, Montreal, Qc., H3A 2T8, Canada 

b. Nuclear and Particle Physics, University of Athens, Athens 15771, Greece 

ABSTRACT 

Recently determined next-to-leading order sets of polarized parton distributions are used 
to study large-pr PP — ► 7 + X at yfs = 38, 100 and 500 GeV. Certain conversion terms, 
necessary to use the above sets, are determined. It is concluded that, to distinguish be- 
tween the above sets, planned RHIC experiments should be successful, in particular at cm. 
photon pseudorapidity ~ 1, and that a proposed HERA-iV fixed target experiment should 
have rather large accuracy at relatively large xt (= 2pr/^/s). 
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I. INTRODUCTION 



Polarized particle reactions have opened a new domain to test perturbative QCD. In this 
domain, in spite of extensive recent data on polarized deep inelastic scattering [[J, the shape 
and the size of the polarized gluon distribution Ag(x) remains an open question. Important 
progress requires experiments on reactions with polarized initial hadrons, dominated by 
subprocesses with initial gluons, like those planned in or proposed in ||. 

Among such reactions, large-pr direct photon production in longitudinally polarized 
proton-proton collisions occupies a prominent place, hence higher (next-to-leading) order 
corrections (HOC) have been determined HH. The results, based on the then existing 
leading order polarized parton distributions, predicted large asymmetries. 

Recently, however, there have been important developments in two respects. First, the 
two- loop longitudinally polarized splitting functions became known ||. Second, making 
use of them as well as of the recent data jlj, sets of new polarized parton distributions have 
been determined. As before, the sets differ essentially in the input Ag(x). In contrast, 
however, to older sets, due to H, the new sets exclude very large values of the integral 
AG = Jq 1 Ag(x) dx. Thus the predictions of [|]|| must be reconsidered, and this is the 
main purpose of the present work. 

Regarding the extension of the Dirac matrix 75 in n = 4 — 2e dimensions, in our work 
j| we adopted an anticommuting 75 (AC) scheme as well as dimensional reduction (see 
below). The two- loop splitting functions of || are determined in a different scheme, and 
one must calculate the terms to be added to the hard scattering cross section of ||] to 
convert to the scheme of || (conversion terms). An additional purpose of this work is to 
determine these terms. 

In Sect. 2 we determine the conversion terms. In Sect. 3, using three new sets of 
polarized parton distributions, we present detailed results and compare them to ||. Sect. 4 
contains our conclusions. 



II. CONVERSION TERMS 



In a given scheme, the n-dimensional split function P^ a (z,e) is written 



ba\ 



(2.1) 



the function P§ a {z) will be termed e-part. In the same scheme, at some factorization scale 
M, let f a /A(x,M) denote parton distribution. In a different scheme, the corresponding 
quantities are distinguished by primes. Then the parton distributions transform according 
to: ITU 



f' a/A (x, M) = f a/A (x, M) + ?M J X ^f a/A (y, M) 
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where /x is a renormalization scale. Polarized parton distributions, Af a /A(x, M), transform 
in a similar manner in terms of the polarized e-parts, AP§ a (z). 

In view of the multitude of schemes, in particular for problems of polarized particles, 
we specify those relevant to the present work. 

The first is the anticommuting 75 scheme previously mentioned. This was introduced 
in and was used to determine the HOC of || and most of the HOC of [|J. In this 
scheme, in addition to AP^(z, e) (= P^(z, e)), we need and have determined AP" s (z,e) 
and AP qg (z,e). The complete expressions are: 



AP? q (z,e) = C F ^^-e(l-z) + ^6{l-z)\ 

AP? g (z,e) = 2N c {y^-2z+l} + (b + s^)6{l-z) 

AP"(z,e) = 2N F (z--) (AC scheme) (2.3) 
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where C F = 4/3 and Nc = 3 (colour SU(3)), N F is the number of flavours and b = 
(lliVc - 2iV»/6. 

The second scheme is dimensional reduction (RD) ||, in which part of the HOC of || 
were determined. In RD, for all a,b: 

P* b (z) = AP^z) = (RD scheme) (2.4) 



The third is the t'Hooft-Veltman (HV) scheme | 1Q| , apart from the e-part of AP g ™, 
which is taken as in ( p.3| ) to satisfy helicity conservation. Here we need all split functions, 
and the remaining ones are found to be (see also the second of 0): 



AP r g ^e) = 2N c \—-—-2z + l + 2e{l-z)\ + {b + e-^5{l 



AP-(z,e) = 2N F \z-±-e{l-z)} 

AP£ q (z,e) = C F {2- z + 2s(l- z)} (HV scheme) (2.5) 

It is this scheme to which we should convert our results. 

The conversion terms are known to arise from differences in the £-parts of the split 
functions and are completely determined from the factorization counterterms 0. The 
form of the latter has been given in ||, and the determination of the conversion terms is 
straightforward. 

Denoting a given subprocess by a(pi) + b(p 2 ) — > j(p) + c + d, where the quantities in 
parentheses are 4-momenta, we define: 

s = ( Pl +p 2 ) 2 i=(pi-p) 2 u=(p 2 -p) 2 
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and 

v = l + i/s w = -u/(s + i) (2.6) 

Below we give the terms to be added to Af = Af(v,w) defined by: {Eq. (4.4) of [4(b)]} 



Ada <3> f V2 f 1 a* 

E— = — dvv(l-v) / dwwAF a/A (x a )AF b/B (x b ){AB(v)S(l-w) + -^Af} 
a 6 p npf, J Vl Jw 1 In 

+ (A <-> B, V <-> -rj) (2.7) 



where Vi,V2,Wi,x a = x a (v,w),x b = x b (v,w) and AB are given in [4(b)]; r] is the cm. 
pseudorapidity of the photon and $ is also specified below. 

First consider g(pi) + q{j>2) lip) + Q + 9- We need the conversion term from AC to 
HV scheme. With $ = iraa s e 2 /N c , writing 

A = v(l-v)w/2, 

and setting 2N F — > 1 in the expressions of AP qg (z,e) we find: 

Af qg = -A-\l - w){C F [v 2 + (1 - v) 2 } + 2N C (1 +v)(l- v) 2 } (2.8) 

For gg — > 7gg we should also go from AC to HV scheme. With $ = 7raa s e^/8: 

Af gg = 2C F A- 1 v 2 (l-w)(3-v + vw) (2.9) 

Now consider <^(pi) + qp{p2) - ► 7(p) + Qa + 9/3- We need the conversion term from RD 
to HV scheme. Here the e-part of the split function P™ q in the (conventional) dimensional 
regularization scheme is required. We find: 

P^ q (z, e) = [1 + (1 - zf]/z - ez (dim. reg.) 

Then, with $ = T\Cpaa s j2N c ^ we find for (3 7^ a: 

„ f 2N 2vw l - w v 3 w 2 (2 - vw), 9r , 2(1 -u) 

A/ OT = 2?;e 1 + ^ + el + ^ 

l — v A (1 — vw) A vw 

1 — w 
vA 

and for j3 = a and z = 1 — v + vw: 



(1 + «)(!-«)=»]} (2.10) 
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Finally we consider q(p\) + q(p2) — > l(p) + Q + q (q, q of identical flavour) and give the 
conversion term from RD — > HV scheme. With $ = 7rCpaa s /2Nc'- 



Af g9 = 2ve 2 {l + 



2 2vw (1 — v) 2 + v 2 w 2 2 v 2 w 2 1 — w ,v 3 w 2 (2 — vw) 
q ^ + 1 - v 7 2 N^(l-v)z (1-vw) 2 



y ,y J -]} (2.12) 



The contributions of qq — > jgg and q a q a — > iqpqp with q^g^ produced via (7 — >■ g^gg 
are obtained from the corresponding unpolarized subprocesses by a change of sign. As in 
[Q, we use unpolarized results in dimensional regularization (MS scheme) and in view of 
taking AP qq = P qq (see before Eq. ( |2.5| )), conversion terms are absent. 



Note that (|2.9|) is the same as the corresponding term of ||11|| , a fact easily understood: 
The terms of jTIJ convert from the so-called MSp scheme, introduced in 0, to HV. The 
MS P scheme has AP e qq = C F [-(1 - z) + 5(1 - z)/2], AP £ gg = N F 5(1 - z)/Q and AP e ab = 
for the rest. Thus, as far as we are concerned, it is equivalent to the AC scheme. Likewise, 
in ( [2.10| )-( |2.12j ), the parts ~ 1/A are the same as the corresponding terms of fllfl ; they 



amount to converting from MSp to HV. Finally, with the transformation v' — 1 — vw, w' 
(l — v)/(l — vw) [i.e. t <-> u], and taking into account that the Jacobian d(v', w')/d(v, w) 



v/(l — vw), we find that (2.9) becomes the same as the corresponding result of [p| ; it 
should be so, since the latter interchanges pi and P2 between initial partons. 

Some doubts have been expressed about the possibility to convert from AC to HV 



scheme [01. In view of the above, and as far as we can tell, the doubts are unfounded. 



III. RESULTS AND DISCUSSION 



We use various sets of polarized parton distributions of one group (the next-to-leading 
order sets A, B, C of fl2]| ); different groups proceed with different input assumptions, and 



this may obscure the degree of real difference in the gluon distribution. 

Furthermore, in presenting asymmetries we always divide the polarized by unpolarized 
cross sections determined via one set of unpolarized parton distributions, of STEQ4M ||13|| ; 



dividing by cross sections determined via distributions of different sets obscures to some 
extend the differences in Ag(x). CTEQ4M has A4 = 0.296 GeV, which we also use in 
our (two- loop) expression of ct s (ix), varying it as we cross the 6-quark threshold (nib = 4.5 
GeV). 

An explanation of our attitude towards 7 Brems and the related 7 fragmentation is in 
order. First, the related effects are known to be important at rather small xp(= 2p^j ' \fs). 
The RHIC experiment, being a colliding beam one, is expected to observe isolated 7's. 
Then, unless isolation criteria are specified, a calculation of 7 Brems is of dubious value. 
The proposed HERA-iV experiment 0, being a fixed-target one, may observe non-isolated 
7's. However, it's energy (yfs = 39 GeV) is relatively low and the values of x F for which 7 
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Brems is important correspond to pt rather small. At small pt, other effects, like intrinsic 
parton's k?, higher twist etc also become important. In particular at lower y/s, the cross 
sections are steeper and kx effects are stronger. 

We thus prefer to leave out 7 Brems; as in [|J, our calculations include only the fac- 
torization counterterms necessary to cancel the mass singularities of collinear 7 emission. 
The subsequent results correspond to xt > 0.08. 

We briefly comment on the input Ag(x) of fL2| (Ql = 4GeV 2 ). In sets A, B, Ag(x) > 
throughout; in set C, Ag(x) changes sign, and for x > 0.1 becomes negative. The integral 
AG has its largest value for A and its smallest for C; even for A, however, it is significantly 



smaller than the large AG of [14] used in [|J. 

As in ||], we present results for pp — > 7 + A at y/s = 38, 100 and 500 GeV; the first value 
is relevant to the HERA- A experiment ||. Also, we consider photon cm. pseudorapidities 
77 = 0, 1 and 1.6 and use /i = M = px- As there are similarities in the pattern, in Figs. 1 
and 2(a), (b), we present results for set A and r\ = 1, but we comment on the results for 
other sets and 77's. 

Beginning with qg — > jq, denote by cr B (gq) [<Jho(.9<i)} the contribution of AB [Af] to 
EAda I 'd 3 p, Eq. ( [2.7|) . Compared to 0], in the considered range of xt, for sets A, B the 
A-factor K gq = [asigq) + ^#0 (#<?)] /^b (<?<?); changes very little; this is clear in Fig. 1(a). 
For set C, at yfs = 38 and 100 GeV and r\ = 1.6 and 1, K gq is not smooth; at r\ — and 
for yfs = 500 at all r], it is similar to Fig. 1(a). 

Next, for gg — > •yqq consider the corresponding Af and denote by a(gg) its contribution 
to EAda/d 3 p. The ratio K gg = cr(gg) / UB^qg), used as a measure, for sets A, B is negative 
and in magnitude very small over all our kinematic range; Fig. 1(b) presents a typical K gg . 
For C it is similar, except for yfs = 38 and 100 GeV and small xt, where K gg > 0, but still 
very small. 

For qq — > jqq, denoting by a(qq) the contribution of the corresponding Af, we consider 
K qq = cr(qq) / UB^qg)- For all sets, K qq is negative. In magnitude, for sets A, B it is very 
small, for C somewhat larger. Fig. 1(c) shows a typical K qq . 

Now we consider the A'-factor for all 0(a s ) and 0(a 2 s ) contributions: 

K = [(T B (gq) + VHo(gq) + cr B (qq) + (?Ho(qq) + cr(gg) + cr(qq)]/[a B (gq) + a B {qq)\ (3.1) 

For sets A, B, due to smallness of \K gg \ and \K qq \, K is almost the same as K qg ; Fig. 2(a) 
makes this clear. As for K qg , for set C, K is smooth only at yfs = 500 GeV. 

Turning to cross section EAda/d 3 p, for sets A, B the shapes are similar to Q but the 
magnitudes smaller; this is clear in Figs. 2(b) and 2(c). We remark [Fig. 2(c)] that, as in 
[|J, the maximum of EAda/d 3 p is at rj > 0, near rj = 1 at the higher energies. For set C 
the behaviour is complicated: at y/s = 38 GeV for 77 = 0, 1 and 1.6, at yfs = 100 for 77 = 
and 1 and at yfs = 500 for 77 = 0, EAda/d 3 p changes sign as xt varies. 

Figs. 3,4 present in detail asymmetries 

A LL (p T ,s,r]) = EAda{p T ,s,rj)/d 3 p/Eda{p T ,s,ri)/d 3 p (3.2) 
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Considering first \fs = 38 GeV, for pt < 6 GeV, due to the smallness of AG, All is small 
for all sets. It becomes larger at larger p?, and for iy > 0.4 there is a clear distinction, at 
least between set A and sets B, C. The cross sections, however, become smaller and the 
proposed HERA-iV experiment || will need high accuracy. At ^Js = 100 and 500 GeV 
(RHIC Spin), to distinguish between A and B, C should not be difficult at either r\ ~ or 
T] ~ 1; to distinguish between B and C, data near rj ~ 1 (Fig. 4) may be necessary. 

Note that, in general, All determined with HOC and without (only Born) differ. Fig. 4 
shows the latter for set A, yfs = 100 and 77 = 1 (dash-dotted line); at the higher pt, it 
differs by a factor ~ 2. 

Finally we turn to the effect of changing the scales, and consider the ratio (Fig. 5, set 

A): 

r = EAda(M = fi = p T /2)/d 3 p/ EAda{M = ft = 2p T )/d 3 p (3.3) 

As one expects, with the evolution of the distributions determined via two-loop splitting 
functions, there is more stability than in With set B, r is even more stable. Set C leads 
to unstable r, ranging from negative to positive values. 

One may wonder about the effect of the conversion terms (Sect. 2) on our results. For all 
quantities considered, the effect does not exceed 6%; of course, mathematical consistency 
requires the presence of the terms. 



IV. CONCLUSIONS 



arc 



Our essential conclusions, based on the next-to-leading order sets A, B, C of [|12 
as follows: 

In general, there is a marked difference between the predictions of sets A, B (Ag(x) > 0) 
and of set C (Ag(x) changing sign). 

For sets A, B, the i^-factors ( |3.1| ) exceed unity; thus the HOC enhance the Born cross 
sections. For set C, fT-factors are more complicated and in several cases the HOC are of 
opposite sign to the Born. 

Due to smaller AG, the cross sections and asymmetries All are smaller than in Q. 
All become significant at relatively large pt, and the differences between All of the above 
three sets are larger at 77 ~ 1, where also polarized cross sections are larger, in general. 
Thus, to distinguish between A, B and C, RHIC, with its expected high luminosity, should 
be successful, in particular near 77 = 1; on the other hand, HERA-iV will need rather high 
statistics at xt > 0.4. 
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FIGURE CAPTIONS 

Fig. 1. Relative contributions of the basic subprocesses using the distributions of set A at 
r] = 1. Dashed lines: yjs = 38 GeV, dotted: y/s = 100, solid: y/s = 500. (a) The 
fT-factor K gq for gq — > jq. (b) The factor K gg = o-(gg)/a B (gq) for gg — > •yqq. (c) The 
factor K qq = a(qq) / 'cr B (gq) for qq — >■ ^qq. 
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Fig. 2. Results with set A: (a) The X-factor ( |3.1| ) for all 0(a s ) and 0(a 2 s ) contributions at 
rj — 1. (b) Inclusive cross sections for pp — > 7 + X vs sj- = Ip^j \fs for 77 = 1. Lines 
as in Fig. 2(a). (c) Inclusive cross sections for pp — > 7 + X vs 77. 

Fig. 3. Asymmetries for all three sets of distributions at pseudorapidity 77 = 0. 

Fig. 4. As Fig. 3 for 77 = 1. 

Fig. 5. Ratio r = [£Ad<x(M = // = p T /2)/d 3 p]/[EAda(M = /j = 2p T )/d 3 p] with set A. 
Solid lines: 77 = 1.6. Dotted: 77 = 1. Dashed: 77 = 0. 
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